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We theoretically design a graphene-based all-organic ferromagnetic semiconductor by terminating 
zigzag graphene nanoribbons (ZGNRs) with organic magnets. A large spin-split gap with 100% spin 
polarized density of states near the Fermi energy is obtained, which is of potential application in spin 
transistors. The interplays among electron, spin and lattice degrees of freedom are studied using 
the first-principles calculations combined with fundamental model analysis. All of the calculations 
consistently demonstrate that although no d electrons existing, the antiferromagnetic n — iv exchange 
together with the strong spin-lattice interactions between organic magnets and ZGNRs make the 
ground state ferromagnetic. The fundamental physics makes it possible to optimally select the 
organic magnets towards practical applications. 

PACS numbers: 73.22.Pr, 75.75.-c, 71.15.Mb 



Introduction. — Since the experimental discovery of 
graphene in 2004 [lj, graphene-based transistors have 
rapidly developed and are considered good candidates 
for post-silicon electronics [2j. However, they encounter 
an outstanding challenge of how to open a sizable and 
well-defined band gap in graphene to control the ON- 
and OFF-states of the transistors. One realistic solu- 
tion is constraining large-area graphene in one dimension 
to form graphene nanoribbons (GNRs) which have band 
aps approximately inversely proportional to their widths 
In recent years, this scheme has been experimentally 
proved rather practical and thus motivated a new frontier 
area called "graphene nanoribbon electronics" [4J46(. 

Besides the gapped semiconductor properties suitable 
for electronics, graphene nanoribbons possess fascinat- 
ing magnetic properties for spintronics. Zigzag graphene 
nanoribbons (ZGNRs) have spatial spin-polarized ground 
states with spin moments coupled ferromagnetically 
(FM) on the same edge and antiferromagnetically (AFM) 
between different edges [3, |8| > which makes them poten- 
tial materials for carbon-based spintronic devices avoid- 
ing heavy magnetic atoms. By means of electric field con- 
trol, Boron/Nitron doping or edge termination, ZGNRs 
can exhibit half-metallicity that may be used for spin 
injection and filtration in metallic spintronics [9r|l2j. 

On the other side, semiconductor spintronics has been 
the subject of many recent studies [13| . The central issue 
is how to develop spin transistors with appropriate ferro- 
magnetic semiconductor materials for non- volatile mem- 
ory applications [14J . By virtue of the unique electronic 
and magnetic properties mentioned above, GNRs may be 
ideal materials for spin transistors. However, the diffi- 
culty is that the ferromagnetic state of ZGNRs is not the 
ground state, and even worse for transistors, it is a gap- 
less state. Therefore, the organic integration of electronic 
and magnetic properties of graphene for spin transistor 
is still an open question. 

In this paper, we design a graphene nanoribbon ter- 



minated with organic magnets to give a possible solu- 
tion for the above issue. In order to convert the ground 
state of ZGNRs from the AFM state to a gapped FM 
one, we use edge passivation of graphene nanoribbons by 
organic ferromagnetic radicals (one kind of organic mag- 
nets, see Fig JTJ). The physics is that the magnetic order 
of the host material will be changed when it exchange- 
couples to the localized spin of attached organic mag- 
nets [15| . Comparing with ZGNRs adsorbing magnetic 
transition metal atoms on their surfaces in literatures 
[la . Il7 |. these carbon-based materials have the advan- 



tages of much weaker spin-orbit and hyperfine interac- 
tions, which should induce much longer spin coherent 
length. Furthermore, owing to the planar structure of the 
selected organic radical, above mentioned exchange cou- 
pling is tunable by manually rotating the radical plane to 
affect the spin polarization near the Fermi energy. Our 
results suggest that the controllable graphene-based spin 
transistors can be expected. 

Methods. — Towards the practical applications, we 
firstly choose trimethylenemethane (TMM) as the or- 
ganic magnet, for the reason that it has been widely in- 
vestigated and designed for kinds of possible spintronic 
applications due to its high magnetic moments [18]. How- 
ever, as will be demonstrated later, the dependence of the 
physics on this specific choice is quite weak. TMM con- 
tains 4 carbon and 6 hydrogen atoms with 7r-conjugated 
structure and has a ground state of triplet biradical, as 
shown in FigJTjD. Two unpaired electrons in TMM pro- 
vide 2fiB spin moments in total distributing via molecu- 
lar orbital over all the carbon atoms. At its spin triplet 
ground states, three outer carbon atoms have the same 
kind of net spin density, while the inner one has the 
opposite (see FigJTJi). When attached to one of the 
ZGNRs' edges already possessing spin moments, TMM 
will strongly modify the properties of ZGNRs. 

6-ZGNR terminated by hydrogen atoms is used as pris- 
tine ZGNR as shown in FigQJa. The unit cell of 6-ZGNR 
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Figure 1: (color online) Top view structures of (a) pristine 
6-ZGNR, (b) TMM, (c) TMM terminated 6-ZGNRs, and (d) 
the net spin density of TMM. The yellow balls are carbon 
atoms and the cyanic ones are hydrogen atoms. The red/blue 
arrow indicates the up/down spin component. 



contains 12 carbon atoms and 2 hydrogen atoms (see the 
figure inside the green rectangle in FigJTJi). The length 
of periodic direction of the unit cell is 2.46 A (experi- 
ment value) , while the repeated images of the molecular 
are separated by vacuum more than 8 A thick along the 
other two directions. To avoid direct interaction between 
images of terminated TMM, the zigzag direction of the 
unit cell is doubled (see the region inside the red rect- 
angle in Fig[T^,). In this supercell, TMM replaces one of 
the four hydrogen atoms, making the system terminated 
by TMM and hydrogen alternately on one edge (e.g. the 
right edge) . The top view of one possible relaxed struc- 
ture (that will be proved the ground state) is shown in 
FigOJ. Please be noted that the TMM biradical is not 
coplanar with ZGNR after atomic optimization. 

The first-principles calculations are performed using 
Density Functional Theory (DFT) method implemented 
with SIESTA code[19| using the Perdew-Burke-Ernzerhof 
exchange-correlation functional [20] . Double-^ plus polar- 
ization function (DZP) basis set is used, together with a 
mesh cutoff of 200 Ry and norm-conserving pseudopo- 
tentials. All the atoms in the supercell are fully relaxed 
to fulfill the energy and force convergence of 10~ 5 eV 
and 0.01 eV/A respectively. The optimized structure 
parameters with different initial spin configurations are 
crosschecked using a plane- wave method implemented in 
Quantum Espresso code|21|. with the same convergence 
criteria and a plane wave kinetic energy cutoff of 30 Ry. 
The identical results ensure the validity of the simula- 
tion parameters used, and other data are handled using 
SIESTA code solely. 

Results. — In what follows, we use "++" and "+-" to re- 



spectively denote the FM and AFM spin configurations, 
where the first symbol denotes the left edge spin moment 
while the second denotes the right one of the graphene 
ribbon. Without losing generality, we designate "+/-" 
to represent up/down spin in what follows. The con- 
vergence of the simulation parameters is firstly checked 
by electronic property calculations of a pristine 6-ZGNR 
with a 1x1x2 supercell. The "+-" order is confirmed to 
be the ground state with a net spin moments of 0.26 /is 
localized on each edge atom, and its total energy is 26 
meV/supercell lower than the "++" one. These values are 
consistent with previously reported results and used as 
reference data [7J, [8j . 

The TMM units are then attached to the right side 
of the ZGNR. Since the spin moments on TMM and 
ZGNRs can align either parallel or anti-parallel, the "++" 
and "+-" spin configurations will extend to four possible 
ones: "+-+", "++-"/'+--" and "+++", with the rightmost 
symbol representing the net spin moments on TMM. Full 
atomic relaxation calculations are done separately using 
these four possible spin configurations to find the ground 
state. It is found that both the bond length and the dihe- 
dral angle between the ZGNR and TMM are varied with 
different spin configurations ( we will discuss this point in 
details later). The relative total energy ( taking E+.+ as 
reference), converged net spin moments and optimized 
bond structures are summarized in Tab|I] Crosscheck- 
ing with Quantum Espresso presents highly consistent 
results, as shown in Tab|I] 

As presented in TabJH the "+-+" state is the ground 
state. It has a total energy of 36 meV lower than that 
of "++-" state, 578 meV than that of "+--" state and 
612 meV than that of "+++" state, respectively. These 
data indicate that, in the ground state of the radical ter- 
minated graphene nanoribbon, the net spin moments of 
TMM radical antiferromagnetically couples to the near- 
est edge spin of ZGNR. As the energy difference between 
pristine "+-""++" ZGNR is 26 meV, the present coupling 
is much stronger than the interedge superexchange of 
ZGNR[22]. A relative result was reported just recently 
by Atodiresei et al. focusing on the spin polarization of 
nonmagnetic benzene, cyclopentadienyl radical and cy- 
clooctatetraene molecule adsorbed onto a ferromagnetic 
2 ML Fe/W(110) surface. It is found that at the organic 
molecule site an inversion of the spin polarization occurs 
with respect to the ferromagnetic surface, resulting from 
the antiferromagnetic coupling between 7r-electrons of the 
molecule and the ci-electrons of Fe atoms [23]. While in 
the present case, the strong antiferromagnetic coupling 
presented here results from the 7r-7r interactions. 

Comparing the corresponding structures of the four 
spin configurations manifests the close connection be- 
tween the spin coupling and the variation of the ZGNR- 
TMM bond length (or dihedral angle). It can be seen 
from TabU that the changes of bond lengths and spin 
moments among the four possible spin configurations 



Table I: Relative total energy, converged net spin moments on selected carbon atoms, optimized bond lengths and spin angles 
for all possible spin configurations of TMM terminated 6ZGNR. The meanings of the symbols for the carbon atoms are depicted 
in FigQJ. The values in brackets are corresponding results obtained by Quantum Espresso code. 
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mainly take place around the C 9 i-Cji bond. When TMM 
and ZGNR are coupled antiferromagnetically ("+-+" and 
"++-" states), the length of the C 9 i-Cti bond is short- 
ened about 0.1A relative to those of "+--" and "+++" 
states. It is reasonable to attribute such a large distortion 
to the softness of the all-organic system and the strong 
spin-lattice interaction. It is further found that, if the 
dihedral angle is disturbed, the bond length will changes 
correspondingly. As the 7r electron orbital trends to be 
perpendicular to the a bonds surface, the dihedral angle 
reflects the spin angle of 7T-7T electrons on C 9 i and Cti. 
Therefore, the calculation reveals that the ir-ir spin cou- 
pling depends not only on the bond length d but also on 
the spin angle 9. We thus propose a spin-lattice (spin- 
phonon) coupling model to illustrate the physics as fol- 
lows, 



H T = H c + J (l - au)S g i ■ Sn + 7« 2 



(1) 



where H c denotes the invariant part (except the spin-spin 
coupling between C 9 i and Cti) before and after spin- 
lattice interaction involved. The second term is the 7r- 
■k spin interaction modulated by the bond distortion u 
between C 9 i and Cti, where Jo is the coupling constant 
at u — 0, a the spin-lattice coupling constant (a > 0). 
The third one is the elastic potential energy caused by 
the bond distortion. As the magnetism of the terminated 
radical is robust, we take a classic approach with 

S gl -S n = S gl Sticos6 (2) 

By taking the differential of the Hamiltonian (fTJ) with 
respect to the bond distortion u, we can obtain the stable 
state under certain fixed angle 9, 

aJoSgiSti cos 9 



Uq cos ( 



(3) 



which indicates that the optimal bond distortion u will 
change with the spin angle 9. We then obtain the optimal 
spin angel 9 and bond distortion u by minimizing the 
total energy, which is given as 



It derives that the ground state is characterized by an- 
tiferromagnetic exchange interaction and a shorter C 9 i- 
Cti bond length, which is consistent with our first- 
principles calculations as shown in Tab. HI 

In order to elucidate the robust physics of the spin- 
lattice coupling, we substitute the TMM radicals with 
CH2, the simplest carbon based organic radical hold- 
ing the sp 2 hybridization. In the first principle calcu- 
lation, the optimal bond length d between C 9 i and Cti 
for each 9 is obtained by minimizing the total energy. 
The calculated d — 9 relation is shown in Fig. [2J which 
is well fitted with a cosine type line as deduced from our 
model. Therefore, in our proposed all-organic ferromag- 
netic graphene structure, the ferromagnetism is provided 
by the n — -k spin interaction, while the spin-lattice cou- 
pling is vital for its stability. 
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Figure 2: (color online) d — 6 relation of the total energy 
minimum state for given spin angle 9. The black squares are 
calculated using DFT, and the blue line is the fitted line using 
equation d=do + u=do + uocos(9, where the fitting parameters 
d =1.470A and u =0.093A. 

Besides the all-organic ferromagnetic ground state 
mentioned above, we find that a large spin-split band 
gap and the 100% spin polarization near the Fermi en- 
ergy are specific characters of TMM and CH2 terminated 
ZGNR. FigJ3^ and b respectively plot their spin-resolved 
band structures and density of states (DOS). From these 
figures, TMM and CH 2 terminated ZGNRs exhibit the 



following similar features: 1) both of the spin up and 
down subbands manifest semiconductor character with 
energy gaps ~ leV; and 2) the spin splitting between 
up and down subbands near the Fermi energy (Ep) is 
about 0.5 eV, which induces a gapped spin-split DOS 
with 100% spin polarization within a wide energy region 
near Ep. The highest occupied and lowest unoccupied 
energy level are proved to be extended states which can 
serve a large current in graphene-based transistors. In 
Figj3j:, the local density of states (LDOS, with energy 
region E-p — > Ep — 0.6 eV) of the highest (spin up) oc- 
cupied energy level of CH2 terminated ZGNR is shown. 
Towards the practical transistors, those features provide 
well-defined conducting ON- and insulating OFF-states, 
and also the high spin-polarized current at the ON-state. 
By applying a positive or negative gate voltage, one can 
selectively shift the (purely spin up) occupied highest en- 
ergy level or the lowest (purely spin down) unoccupied 
one towards the Fermi energy. A 100% spin polarized 
current is thus produced via the gate voltage control. 

For comparison, we make some comments on the 
nonmagnetic side group terminated ZGNRs in the 
liter ature[l2j. For those materials, the sizable energy 
gap and 100% spin polarization in large energy range 
near Ep can not obtained simultaneously. For example, 
in the NH 2 , OH, and COOH terminated ZGNRs, the en- 
ergy gap is about 0.3 eV, but the spin splitting of the 
states near the Ep is very small. In the NO2 terminated 
ZGNR, although the spin splitting is relatively enlarged, 
the energy gap reduced to about 0.1 eV, too small to be 
used in transistors. Therefore, the present design has an 
advantage in realizing a large energy gap and a high spin 
polarization near the Fermi energy simultaneously. 

In spintronics, the spin polarization P{E) may be de- 
fined as 



P(E) = 



DOSf(E) - DOSi(E) 
DOS t (E)+DOS l {E)' 



(5) 



which can be tuned by changing the spin polarized DOS 
near Ep. For this purpose, we study the evolution of the 
spin polarization near Ep with manually varying the spin 
angle 9 in CH2 terminated ZGNR. The result is demon- 
strated in Fig. [5Ji. As mentioned above, the DOS above 
and below the Ep are both 100% spin polarized in a wide 
energy range at 9 — it. However, as 6 decreases towards 
7r/2, the high spin polarized region near Ep shrinks and 
even almost disappears at 9 = tt/2. As seen in the figure, 
the large value of DOS near Ep during radical rotation 
can keep the ON-state under certain gate voltage. Thus, 
controlling the exchange coupling strength by radical ro- 
tation is an effective way to tune the spin polarization. 

Finally, let us comment on the possibility of taking 
Lieb's theorem on the repulsive Hubbard model [24| as 
the mechanism of the gapped ferromagnetic state, for the 
reason that its deduced Klein's edge [25j has been applied 
to explain the ferromagnetic state of some other edge ter- 
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Figure 3: (color online) (a)(b) Spin-resolved band structure 
(left panel) and density of states (right panel) of the ground 
state of TMM and CH 2 terminated ZGNR, respectively, (c) 
Local density of states(LDOS) of the ground state of CH2 
terminated ZGNR. Red for up spin, (d) Spin-resolved DOS 
of CH2 terminated ZGNRs with spin angle 8 equal to tt, 37r/4, 
27r/3, and it/2 respectively. The shaded area presents the spin 
polarization P(E) defined in Eq. ([5]). 



minated GNRs in literatures [7J, 126( . The central point of 
the Lieb's theorem is that the ground state of a bipartite 
lattice and half-filled band has a net total spin S provid- 
ing the two sublattice have different number of sites. By 
carefully comparing our results with Lieb's theorem, we 
conclude that it can not consistently explain all the fea- 
tures of the ground state here. The main discrepancies 
are summarized as follows: 1) it is the strong spin-lattice 
interaction, rather than the Hubbard (electron-electron) 
correlation, stabilizes the ground state of the present sys- 
tem; 2) the spatial spin order of the ground state is be- 
yond Lieb's theorem that can only deduce a nonzero total 
spin; and 3) our DFT and model calculations have both 
omitted the quantum fluctuation of electron-electron in- 
teraction which is the curial part of Hubbard model and 
Lieb's theorem. Relative to Lieb's theorem, our mech- 



anism of the antiferromagnetic exchange together with 
strong spin-lattice interactions can explain all of the fea- 
tures of the gapped ferromagnetic state in a consistent 
manner, as already demonstrated above. 

Summary. — In summary, we designed a graphene- 
based all-organic ferromagnetic semiconductor by termi- 
nating ZGNR with organic mangets. A large spin-split 
gap with 100% spin polarized DOS near E-p is obtained, 
which is of potential usage in spintronic devices. Combin- 
ing first-principles calculations with basic model analysis, 
we conclude that the mechanism for the gapped ferro- 
magnetic state is the antiferromagnetic 7r — 7r exchange 
together with strong spin-lattice interactions between or- 
ganic magnet and graphene nanoribbons. By controlling 
the 7r — 7r interaction strength, the spin polarization of 
the DOS near Ep can be tuned from 100% to nearly zero. 
The fundamental physics makes it possible to optimally 
select the organic magnet towards practical applications. 
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